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^^ ' Aims. We present the results of contemporaneous spectroscopic and photometric monitoring of the young solar-type star 
r^ ^ HD 171488 (Prot — 1.337 days) aimed at studying surface inhomogeneities at both photospheric and chromospheric levels. 

Methods. Echelle FOCES spectra {R ~ 40 000) and Johnson BV photometry have been performed in August 2006, with a 
good coverage of rotational phases. Spectral type, rotational velocity, metallicity, and gravity were determined using a code 
developed by us (Rotfit) and a library of spectra at the same resolution of slowly-rotating reference stars. The metallicity 

Q ■ was measured from the analysis of iron lines with the MOOG code. The spectral subtraction technique was applied to the most 

}~i ' relevant chromospheric diagnostics included in the FOCES spectral range, namely Call IRT, Ha, HeiDa, H/3, and Can H&K 

C/3 , lines. 
, ^, ' Results. A simple model with two large high-latitude spots is sufficient to reproduce the B and V light curves as well the radial 
velocity modulation, if a temperature difference between photosphere and spots of about 1500 K is used. A Doppler imaging 
^—( , analysis of photospheric lines basically confirms a similar spot distribution. With the help of an analogous geometric two-spot 

^ ■ model, we are able to reproduce the observed modulations in the residual chromospheric emissions adopting different values of 
^^ I ratios between the flux of plages and quiet chromosphere (about 5 for Ha and 3 for Call diagnostics). Facular regions of solar 
^^D type appear to be the main responsible for the modulations of chromospheric diagnostics. Both the spot/plage model and the 
'•'i cross-correlation between the light curve and the chromospheric line fluxes display a sensible lead effect of plages with respect 
^— ^ to spots (from 20° to 40° in longitude), as already observed in some active solar-type stars and RS CVn systems. 
l/^ • Conclusions. The contemporaneous monitoring of photospheric and chromospheric diagnostics in the young and rapidly rotating 
c^J ' solar-type star HD 171488 allowed us to detect active regions which have nearly the same location at both atmospheric layers, 
^^ , with plages slightly leading spots in longitudes. These active regions are similar to the solar ones in some respect, because the 
^~~^ • spot temperature is close to that of sunspot umbrae and the plage flux-contrast is consistent with the average solar values. The 

►^ ' main differences with respect to the Sun are larger sizes and higher latitudes. 
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braking. It has been proposed that the rate of angular mo- 
mentum loss decreases or, at least, saturates at high ro- 
tatio n rates (JMacGreeor fc Brenneiill99ll : iBarnes fc Sofia 
Il996l) . Saturation of magnetic activity is witnessed by the 
typical behaviour of X-ray emission and ot her activity sig- 
natures as a function of the rotation speed (|Randichlll997 ; 
iKrishnamurthi et al.lll997 ). Periodic variability of satu- 
rated fast-rot ating stars is present both in the o ptical con- 
tinuum fe.g-. lMarilh et al.lll997t IStassun et al.lf2004il and 
in the X-ray band (JFlaccomio et al.ll2Q05l ). This rises ques- 
tions about the true extent of saturation in fast-rotating 
stars. It is intriguing that very fast rotating stars, for 
which a "super saturation" regime, i.e. a slight decline 
of magnet ic activity with the increase of rotation rate, is 
displayed (jProsser et al.lll996l ). still show organized mag- 
netic topologies but inefficient angular momentum losses 
(|St(;;pieh et al.ll2004l) . 

For all these reasons we consider it mandatory to inves- 
tigate the surface topology of magnetic active regions on 
fast-rotating young stars, also because of the implications 
on the stability and early evolution of planetary systems. 

HD 171488 (=V889 Her) is a young fast-rotating so- 
lar analogue that has been the subject of several stud- 
ies in the last decade. Its kinematic properties indicate 
that HD 171488 is a member of the Local Association, 
a stream of young stars wi th ages ranging from 20 to 
150 Myr (JMontes et al.ll200ll ). It is probably the brightest 
single early-G type star that is rotating fast enough for 
mapping its photospher e through Doppler Imaging (DI). 



Strassmeier et al.l (|2003f) made the first dedicated photo- 
metric and spectroscopic study of HD 171488 in which 
they defined its astrophysical parameters, concluding that 
it is a single GOV main-sequence star with an age of about 
30-50 Myr based on both its position on the HR diagram 
and the lithium content. They also made the first recon- 
struction of its photosphere through the DI technique, 
showing a big polar spot with additional high-latitude fea- 
tures. Subs equent works based on DI have a lways found a 
polar spot (|Marsden et al.ll200d Ijeffers fc Donati 2008). 
From long-term photometry and DI, two active longi- 
tudes separate d by about 180° seein to be pers istent for 2.2. Photometry 



to the Can H&K (A =3968.49 A, 3933.68 A), He (A = 
3970.074 A), Ha (A =6562.849 A), and Can infrared 
triplet (IRT; A =8498.06 A, 8542.14 A, 8662.17 A) lines. 
The work is organized as follows. In Section[2]we de- 
scribe the observations and data reduction. The spectral 
classification, rotational and radial velocities, evolutionary 
status, and metallicity are discussed in Section|ni The di- 
agnostics of photospheric and chromospheric activity are 
analyzed in Sections |4] and [5l respectively. A short discus- 
sion on the behaviour of photospheric and chromospheric 
activity is presented in SectionlHl Section [7] contains our 
conclusions. 

2. Observations and reduction 

2.1. Spectroscopy 

We observed HD 171488 at the 2.2-m Cassegrain tele- 
scope of the Calar Alto Observatory (CAHA, Sierra de 
Los Filabres, Spain) with th e Fiber Optics Cas segrain 
Echelle Spectrograph rFOCES: IPfeiffer et al.lll998h during 
four nights from 13 to 16 August, 2006. The 2048x2048 
CCD detector Site#ld (pixel size = 24/Ltm) allowed us to 
achieve, with an exposure time of 15 minutes, a signal- 
to-noise (S/N) ratio in the range 80-160 for HD 171488 
{V ~ 7™. 4), depending on airmass and sky conditions. 

The spectral resolution, as evaluated from the full 
width at half maximum (FWHM) of the emission lines 
of the Th-Ar calibration lamp, was in the range 0.15- 
0.22 A from the blue to the red, yielding a resolving power 
i? = A/AA -40000. 

The data reduction was performed by using the 
ECHELLE task o f the IRAFlj packag e, following the scheme 
summarized bv lBiazzo et al.l ( 20091 ). 

To remove the telluric water vapor lines at the Ha and 
NaiD2 wave l ength s, we used the procedure described by 
iFrasca et al.l ( 20001 ). adopting as telluric template a high- 
S/N spectrum of a Cep (A7 IV, wsini = 196 kms~^), 
acquired during our observing run. 



several years (iJarvinen et al.l 120081 : iHuber et al. 2009). 



Jeffers fc Donatil (|2008l ) found a high solar-type differ- 



ential rotati on with the e quator lapping the poles every 
12-13 days. Ijarvinen et al. (2008) found instead a much 
weaker differential rotation. 

Despite the number of works devoted to this object 
in the last few years, which are mainly based on DI or 
Zeeman-DI, a simultaneous detailed investigation of the 
chromospheric and photospheric inhomogeneities has not 
been performed so far. 

In the present work we analyze contemporaneous BV 
photometry and high-resolution spectra of HD 171488. We 
determined the main parameters of the photospheric spots 
by means of the DI technique applied to a few suitable 
absorption lines and of a spot model for simultaneous so- 
lution of light and radial velocity curves. The behaviour 
of chromospheric inhomogeneities was investigated thanks 



The photometric observations were performed in the B 
and V Johnson filters with the 91-cm Cassegrain telescope 
at the M. G. Fracastoro station [Serra La Nave, Mt. Etna, 
Italy) of the Osservatorio Astrofisico di Catania (OACt). 
The observations were made with a photon-counting re- 
frigerated photometer equipped with an EMI 9893QA/350 
photomultiplier, cooled to —15 °C. The dark current noise 
of the detector, operated at this temperature, is about 1 
count/sec. 

HD 171488 was observed from 14 to 21 August 2006 
for a total of 8 ni ghts, a long with HD 171286 (V = 6'".84, 
B — V= 1.06. IOialll987[ ) used as a comparison star (C) for 



^ IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of the 
Universities for Research in Astronomy, inc. (AURA) under 
cooperative agreement with the National Science Foundation. 
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the differential photometry and HD 171623 {V = 5™. 79, 
B-V = 0.00, present work) and BD+17 3634 (V = 8'°.43, 
B~V — 0.57, present work) as check (Cfci, Ck2) stars. We 
adopted HD 171286 as the local standard for the deter- 
mination of the photometric instrumental "zero points" . 
Several standard stars, selected from the list of Landolt 



(|l992r) , were also observed during the run in order to deter- 
mine the transformation coefficients to the Johnson stan- 
dard system. The observed magnitudes were corrected for 
atmospheric extinction using the seasonal average coeffi- 
cients for the Serra La Nave Observatory. 

The data were reduced by means of the photomet- 
ric data reduction package phot de signed for the pho- 
toele ctric photometry of the OACt (jLo Presti fc Marilh 
1993[ ). The photometric errors, estimated from measure- 
ments of standard stars with a brightness comparable 
to the program stars, are typically cry ~ 0"\010 and 
as-v ~ 0'°.014. 

In order to improve the photometric precision we av- 
eraged four consecutive measurements of the variable and 
adopted the standard deviation, ranging from about 0.005 
to 0.015 mag, as an error estimate for each of the mean 
photometric points. 

The photometric data are reported in Table laFl. 

3. Astrophysical parameters 

Thanks to the high resolution and the wide wavelength 
coverage of the FOCES spectra, we have redetermined 
spectral type, effective temperature, gravity, metallicity, 
radial and rotational velocity and quoted these values in 
Table[T]along with literature values. The evolutionary sta- 
tus of HD 171488 has also been checked. 

3.1. Spectral type and rotational velocity 

The spectral type and the v sin i o f HD 171488 was d erived 
through the IDlJfl code Rotfit (JFrasca et al.ll2006l ). The 
code simultaneously finds the spectral type and the v sin i 
of the target searching for, into a library of standard star 
spectra, the spectrum which is best fitting (minimum of 
the residuals) the target one, after the rotational broaden- 
ing by convolution with a rotational profile of increasing 
V sin i at steps of . 5 km s^ ^ . We acquired spectra of a very 
small sample of standard stars (9 objects) with FOCES 
during our run. Thus, we preferred to use a library of 
185 ELODIE Archive standard stars well distributed in 
effective temperature, spectra l type, and gravity , and in a 
suitable range of metalicities (iPrugniel fc Soubir an 2001) 
that have nearly the same resolution {R ~ 42 000) as our 
FOCES spectra. We found for HD 171448 a G2V spec- 
tral type, a nearly solar metallicity, [Fe/H]=+0.01±0.09, 
log5 = 4.30 ±0.15, and Teff = 5750±130K (see Table [J), 
as average parameters of the best ten standard stars per 



Table 1. Physical parameters of HD 171488 from the lit- 
erature and present work. 



Parameter 



Value 



Reference 



Spectral Type 


GOV 




G2V 




G2V 


V sini (krns^ ) 


38 




39.0±0.5 




37.5±0.5 




38.0±0.5 




37.1±1.0 


Prot (days) 


1.3371±0.0002 




1.313±0.004 




1.33697 


Wlu (mA) 


220 




213±7 




215±7 


Tcff (K) 


5830 ± 50 




5808 




5750 ± 130 




5850±100 


logs 


4.50 




4.30 ±0.15 


[Fe/H] 


-0.5 




-h0.01±0.09 




+0.15±0.07 



Harlan (1969) 
CutisDoto et al. ('2002'! 
Present work 
Fekel fl997) 

^trassmeier_et_^L (2003) 
Marsden et al. (2006') 
Jeffers fc Donati (2008) 
Present work 
^trassmeier_et^L (2003) 
Mars den et al. (2006) 
Jarvinei^_et_^L (2008) 
Cutispoto et al. (2002) 
^trassmeier_et_^L (2003) 
Present work 
Str assmeier et al. (2003) 
Jeffers fc Donati (2008) 
Present work (Rotfit) 
Present work (Moog) 
^trassmeier_et_^L (2003) 
Present work 
Strassmeier et al. (2003) 
Present work (Rotfit) 
Present work (Moog) 



■^ Available as on-line material. 
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IDL (Interactive Data Language) is a registered trademark 
of ITT Visual Information Solutions. 



each echelle order. This temperature range encompasses 
all the previous Toff values from the literature. 

We remark that our temperature determination is not 
strongly affected by starspots, due to their low contribu- 
tion to the observed flux in t he spectral ra, r ige of the obser- 
vations. Indeed, as shown bv lFrasca et al.l (J2005I) . for large 
spot-photosphere temperature differentials the star's inte- 
grated emission is dominate d by non-spot t ed ph otospheric 
flux. In their Appendix A, iFrasca et al.l ( 20051) outline a 
formalism for evaluating the mean temperature observed 
from a star in a given spectral region with a given spot fill- 
ing factor and photosphere-spot temperature difference, 
AT. Adopting T^s = 5800 K, AT = 1500 K, and a fllling 
factor of 0.1, that is likely an overestimation, the measured 
temperature would be, at most, only 50 K cooler than the 
true photospheric temperature. Thus we conclude that the 
"undisturbed" photospheric temperature must be of about 
5800 K and we adopted this value for the following analy- 
sis. 

For the v sin i determination, we chose as "non- 
rotating" templates three stars with spectral type simi- 
lar to HD 171488 observed with FOCES during the same 
run, namely 54 Aql (F8V), 10 Tau (F9V-IV), and 72 Her 
(G2V). We found, as the average of different echelle or- 
ders and these three templates, vsini = 37.1±1.0 kms^^, 
which is in very good agreement with the literature values 
reported in Tablc[T] However, using the Rotfit code with 
the ELODIE templates, we found a nearly equal value of 
37.4kms~i. 
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Examples of the application of the Rotfit code to 
two different spectral regions are shown in Fig. [3q, where 
the very good agreement between observed and standard 
spectra is apparent. 

3.2. Radial Velocity 

We measured the heliocentric radial velocity (Vj-) by 



means of the cross-cor relation technique, (e.g., ISimkin 
1974 iGunn et al.l Il996l) , taking advantage of the wide 



spectral coverage offered by FOCES. The radial velocity 
standard star a Ari (T4 = —14.2 kms~^), observed in the 
same run, was used as template. We cross-correlated each 
spectral order of the FOCES spectra of HD 171488 with 
the template by using the IRAF task FXCOR, avoiding 
the orders with low S/N ratio or contaminated by broad 
and/or chromospheric lines (e.g.. Ha, NaiiD2, Can H&K) 
or by prominent telluric features. We ended up with 60 
orders useful for the calculation of the cross-correlation 
functions (CCFs). The radial velocities listed in Table [5] 
per each spectrum are weighted averages of the values of 
the 60 chosen orders with weights Wi = l/cTi, where ai is 
the error for the i-th order evaluated by FXCOR. The aver- 
age value over the entire observing run is K- — —23.3 ±0.2 
km s~^ , which is in close agreement with previous determi- 
nations. Anyway, the individual values of V^ display a clear 
rotational modula tion (see Fig. E]) as already found by 
Huber et al.l (2009|) and explained in terms of a "Rossiter- 
McLaughlin" effect caused by stars pots instead of eclipses, 
as in the original formulation fe.g.. lRossiteiill924| ). 



3.3. Evolutionary status, lithium abundance and age 

We checked the evolutionary status of HD 171488 by 
means of its position onto the Hertzsprung-Russell (HR) 
diagram (Fig. [T|) and the lithium content. We used the 
evolutionary tracks calculated by IPalla fc Stahlen (|l999r) 
for the pre-main sequence phase. The unspotted photo- 
spheric temperature of 5800 K was adopted. Firstly, to 
derive the stellar luminosity, we evaluated the interstel- 
lar extinction Ay ~ 0.03 mag (quite negli gible) from the 
Hipp arcos parallax (tt — 26.87±0.89mas. lPerrvman et al.l 
19971 ). assuming a mean extinction of 0.7 mag/kpc. Then, 
the de-reddened Vq magnitude at maximum brightness 
(spottedness minimum) was converted into absolute mag- 
nitude My with the parallax and subsequently converted 
into bolometric mag nitude by usin g the bolometric cor- 
rection tabulated by iFloweiJ (|l996h as a function of the 
effective temperatu re. The bo lometric magnitude of the 
Sun, Mboi = 4.64 (jCrog |2000|) , was used to express the 
stellar luminosity in solar units {L ~ I.2L0). 

The compa rison with the e volutio nary tracks and 
isochrones of IPalla fc: Stahlerl (1999) indicates that 
HD 171488 is a post-T Tau or ZAMS star with an age 
of sa50 Myr, a mass of about 1.08 Mq, and a radius 
i? ~ 1.1 Rq , in good agreement with previous determi- 




3.75 

log T^,, 



Fig. 1. Position of HD 171488 on the H R diagram. The 



Palla fc Stabler 



evolut ionary tracks and isochrones of 

(|l999l ) are shown by c ontinuous and dotted l ines, respec- 
tively. The birth line of lPalla fc Stahleij (|l999l) is displayed 
by a dashed line. 
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nations. The comparison with the ISiess et al.l ( 20001 ) evo- 
lutionary tracks leads to slightly larger values of mass 
(M ~ 1.16 Mq) and age (70-100 Myr) that are, however, 
still consistent with the former within the errors. 

A further evidence of the youth of our object is given 
by the presence of the Li i A6707.8 absorption line. Indeed, 
lithium is strongly depleted from the stellar atmospheres 
of late-type stars when mixing mechanisms pull it deeply 
in their convective layers and a deep Li i A6707.8 photo- 
spheric absorp tion line is generally co nsidered as a youth 
indicator (e.g.. lSoderblom et al.lll998l ). 

The equivalent width of the lithium A6707.8 line, 
Wui = 224 ± 8 mA, is the average value of all the EWs 
measured in each FOCES spectrum. We derived a lithium 
abundance Au = 3.2 by interpolation of the NLTE cu rves 
of growth tabulated by IPavlenko fc MagazziJ ( 19961 ). A 
contribution of 9 mA due to the Fei A 6707.4 A line, 
eva luated according t o the empirical correction proposed 
by ISoderblom eraD (|l993l) . A^l^Li(mA) = 20 {B - 
V)o — 3, was subtracted before calculating the abun- 
dance. Although the lithium content cannot be simply 
converted into star age, HD 171488 seems to be slightly 
more lithium-rich than the stars in Pleiades upper en- 
velope of the same temperature. The lithium abundance 
indicates an age in the range 50-150 Myr, in agreement 
with the position on the HR diagram. 

3.4. Metallicity 

Iron abundance was measured, in the local thermody- 
namic equilibri um (LTE) ass umption, using the 2002 ver- 
sion of MOOG (ISnedenlll973h and a grid of ID model at- 
mospheres bv iKuruca (1993). LTE iron abundance was 
derived by means of equivalent widths of 31 Fei and Fen 
lines in the 5000-6800 A range measured using a Gaussian 
fitting procedure with the IRAF task Splot. For the anal- 
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ysis, w e used the prescriptions given by iRandich et al 
( 20061) . We refer to that paper for a wide description of 
the procedure, hue hst, and cr-chpping criteria. 

The effective temperature was determined by impos- 
ing that the iron abundance does not depend on the ex- 
citation potentials of the lines. The microturbulence ve- 
locity ^ was determined by imposing that the iron abun- 
dance is independent of the equivalent widths of Fe i lines. 
The surface gravity log g was determined by imposing the 
Fel/Fell ionization equilibrium. The initial value for the 
effective temperature was the one we adopted as undis- 
turbed photospheric T^s = 5800 K (Sect. |3TT])- For the 
surface gravity, we adopted as starting value the one de- 
rived by means of the Rotfit code for spectral synthesis 
(log (7 = 4.3; Sect. 13. ip . The initial microturbulence ve- 
locity was set to be 1.5 kms^^. Final astrophysical pa- 
rameters were Teg = 5850 ± 100 K, Xogg = 4.3 ± 0.2, 
^ = 1.8 ± 0.2kms~^, i.e. the effective temperature and 
gravity did not change appreciably, while the iron abun- 
dance was [Fe/H]=-f 0.15±0.07 which is in agreement, tak- 
ing the errors into account, with the value provided by 
Rotfit, but it is noticeably higher th an the value of —0.5 
suggested bv IStrassmeier et al.l (J2003l) . 



4. Photospheric activity 

4.1. Light curve and spot modeling 

The photometric data acquired contemporaneously to the 
spectroscopic ones allowed us to obtain a light-curve show- 
ing a rotational modulation due to the presence of spots 
(Fig. [HI . The rotat i onal p hases were derived according to 
the : Jarvinen et al.l ( 20081 ) ephemeris: 



HJD 



0=0 



2 449 950.550-H 1.^^33697 X ^, (1) 



where the rotational period (Pro t) is nearly identical to th e 
value of 1.3371 days found bv IStrassmeier et al.l (|2003[ ). 
With a periodo gram analysis and a CLEAN deconvolu- 
tion algorithm ([Roberts et al.l 119871) we find a period of 
1.37 lb 0.12 days that is fully consistent with the value of 
Jarvinen et al.l (|2008l ) taking the limited time span of our 
data into account. 

The V light-curve has a slightly asymmetric shape with 
a maximum brightness Knax — 7.™40 at about phase 0.4, 
a minimum around 0.0 phase, and a variation amplitude 
Ay :^0.07 mag . Its shape is very similar to that shown by 
Jarvinen et al.l ([2008) in their Fig. 2 for the mean epoch 
2006.79. 

For a first reconstruction of the photospheric inhomo- 

a spot 



geneities, we used Macula (jFrasca et al 



200i), 



model code which assumes circular dark or bright spots 
on the surface of a spherical limb-darkened star. The flux 
contrast between spotted areas and quiet photosphere 
(Fsu/Fuh) is evaluated through the PHOENIX NextGen 
(JHauschildt et al.lll999t ) atmosphere models adopting T — 
5800 K for the quiet photosphere. With only two pho- 
tometric bands and the rather small amplitude of the 



light curve it is not easy to derive the spot temper- 
ature. Thus we preferred to fix the spot temperature 
AT = 1500 K, close to the values f oun d for high lati- 
tude s pots bv lStrassmeier et al.l (|2003l ) and ljarvinen et al. 



(|2008l ) with DI technique, and let the spot area and loca- 
tion be free to vary. In order to minimize the degrees of 
freedom and still allowing the code to reproduce an asym- 
metric light curve, only two spots have been put on the 
photosphere. We a dopted the same inclinat ion of the ro- 
tation axis used bv IStrassmeier et al.l (|2003l ). i = 55°. We 
also exploited the information contained in the radial ve- 
locity searching for a solution of the B, V light curves 
that simultaneously reproduce the VJ- behaviour. The syn- 
thetic curves are displayed by dotted lines superimposed 
to the data in Fig.|6|and the spot parameters are quoted in 
Table m The light curves alone mainly constrains the spot 
longitudes and areas, with only a raw indication of lati- 
tudes. Anyway, the broad light modulation, as well as the 
shape and amplitude of radial velocity variations, force 
the spots to have high latitudes (75° and 60°) as found 
in previous studies. This simple model will be also ap- 
plied to the chromospheric line-fluxes curves in Section [S] 
for the study of the spatial relation between photospheric 
and chromospheric active regions (ARs). 

4.2. Doppler Imaging 

For a more accurate surface reconstr uction we u s ed th e 
Doppler imaging code TempMap bv iRice et al.l ( 19891 ). 
It performs a full LTE spectrum synthesis by solving the 
equation of tr ansfer through a set of model atmospheres 
(|Kuruczl I1993I ) at all aspect angles and for a given set 
of chemical abundances. Simultaneous inversions of the 
spectral lines as well as of two photometric band-passes 
(Johnson B and V in the present case) are then carried 
out using maximum-entropy regularization. 

The 10 available spectra collected over 3 rotations 
cover the whole rotational phase, i.e., they can be used to 
reconstruct one single Doppler image. The last spectrum 
from the time series was excluded due to an insufficient 
S/N ratio. 

Doppler imaging was performed within the 6392- 
6440 A spectral range using two lines, the Fe 1-6411, and 
Ca 1-6439 lines. The quality of the data with a low S/N 
ratio did not allow us to use other lines in the domain. The 
individual reconstructions, as well as the average map, are 
shown in Fig. [5] and Fig. [31 respectively. The maps re- 
vealed a similar spot distribution, i.e., mainly cool polar 
spots with temperature contrasts of up to ssl500K with 
respect to the undisturbed surface of 5800 K. Some low- 
latitude features are also recovered, however, with signif- 
icantly weaker contrast ranging from «300K (Ca 1-6439) 
to a maximum of »i500K (Fe 1-6411). Numerous bright 
features are also appeared, as mirroring of dominant cool 
spots, thus most of them are believed to be artifacts. 

Despite small differences, the resulting Ca and Fe maps 
are in very good agreement (cf. the average map in Fig.|3|). 
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Fig. 2. Doppler maps (top panels), line-profile fits {middle panels), and light-curve fits {bottom panels) for Cai-A6439 
{left panels) and Fel-A6411 {right panels). The arrows below each map indicate the phases of spectroscopic observations. 



The overall structure reveals a dominant cool polar spot 
with a broad appendage extending in longitude between 
~ 180°-300° and another cool region, at a mid-latitude of 
^ 60°, centered at a longitude ~ 45°. This result is rem - 
iniscent of the maps obtained by Ijarvinen et al.l ( 20081 ) , 
where a similar temperature distribution was recovered 
over several years. Moreover, the polar spot/appendage is 
consistent with the high-latitude circular spot (# 2) found 
by modelling our light curve with Macula while the # 1 
region found by the light-curve model at = 10° (see 
Table H]) matches well the lump of the polar spot at ~ 0° 
phase along with the mid-latitude spot at = 45°. 



Ca + Fe combined 
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5. Chromospheric activity 

The wide wavelength range of FOCES allowed us to study 
the chromosphere of HD 171488 by using several lines from 
the near UV to the NIR wavelengths (namely, Ca ii H & K, 
H/3, HeiDa, Ha, Can IRT), which provide information 
on different atmospheric levels, from the region of tem- 
perature minimum to the upper chromosphere. To derive 
the chromospheric losses, we used the "spectral synthe- 
sis" technique, based on the comparison between the tar- 
get spectrum and an observed spectrum of a non-active 
standard star (called "reference spectrum"). The differ- 
ence between the observed and the reference spectrum 
provides, as residual, the net chromospheric line emission, 
which can be integrated to produce an emission equiv- 
alent width, W"^ (see, e. g.. iHerbiel Il985l: lBarde"nlll985 : 



Frasca fc Catalanol 1994; Montes et al.lll995l )~ 



The non-active star used as a reference for the spec- 
tral subtraction is 72 Her (=HD 157214), a G2V star 
{B — V = 0.62) with a very low activity level, as in- 
dicated by the low value of th e Can Mt. Wilson index 
{S = 0.156. lDuncan et al. 1991), close to the value for the 
quiet Sun {S = 0.159, Oranje 1983), and by the deepest 
Ha co re among the G-type stars investigated by iHerbig 
(|l985l ). This star was also observed with FOCES during 
the same run as HD 171488. In Fig. 21 we show an example 
of a spectrum of HD 171488 in the Ha, H/3, HeiDa, and 
Can IRT regions, together with the standard-star spec- 
trum rotationally broadened to wsini — 37 kms~^ which 
mimics the active star in absence of chromospheric activ- 
ity. The Can H & K region is displayed in Fig. [SJ The Ha, 
H/3, and Can IRT profiles are clearly filled-in by emis- 
sion. The Can IRT A8542 line displays a small emission 
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Table 2. Radial velocities and parameters of the subtraeted spectra. 



HJD 

(+2 400 000) 


Phase 


(kms-i) 


(A) 


(A) 


(A) 


TT/cm 
^^Call-IRT 

(A) 


W^C™l-H+K 

(A) 


53 961.346 


0.914 


-23.16+0.58 


0.528+0.061 


0.136+0.026 


0.044+0.020 


1.109+0.073 


0.733+0.076 


53 961.439 


0.984 


-23.37+0.57 


0.575+0.066 


0.141+0.024 


0.054+0.021 


1.140+0.057 


0.764+0.077 


53 961.533 


0.055 


-23.49+0.55 


0.562+0.058 


0.161+0.022 


0.067+0.028 


1.165+0.064 


0.779+0.111 


53 962.355 


0.670 


-23.33+0.57 


0.536+0.066 


0.142+0.022 


0.060+0.023 


1.144+0.051 


0.776+0.081 


53 962.513 


0.787 


-23.49+0.49 


0.469+0.036 


0.118+0.015 


0.062+0.012 


1.065+0.040 


0.722+0.054 


53 962.589 


0.845 


-23.21+0.53 


0.480+0.041 


0.152+0.032 


0.064+0.020 


1.076+0.053 




53 963.465 


0.500 


-22.53+0.57 


0.488+0.042 


0.131+0.018 


0.059+0.026 


1.009+0.058 


0.697+0.054 


53 963.577 


0.583 


-23.00+0.60 


0.484+0.086 


0.164+0.027 


0.085+0.032 


1.032+0.095 


0.709+0.109 


53 964.397 


0.197 


-23.67+0.56 


0.593+0.082 


0.144+0.024 


0.052+0.027 


1.123+0.077 


0.765+0.070 


53 964.536 


0.301 


-23.46+0.53 


0.550+0.061 


0.112+0.030 


0.078+0.023 


1.095+0.085 


0.706+0.084 


53 964.612 


0.358 


-23.17+0.61 


0.511+0.109 


0.149+0.039 


0.091+0.053 


1.114+0.143 





Table 3. Radiative chromospheric losses, y^, and proba- 
bility for the line-flux modulations. 



Line 



Flux 



Fourier pol. 



xl 



X-'<Xo 

Constant 



Hq 

H/3 

He I Da 

CaiiH+K 

CaiilRT 



3.8x10'' 
1.4x10^ 

6.5x10'^ 
5.3x10'' 



0.32 
2.22 
1.53 
0.02 
0.27 



0.3% 6.38 

10% 5.20 

4% 2.77 

0.05% 1.21 

0.2% 5.18 



30% 

12% 

4% 

1% 

18% 



age W^"^ by the continuum surface flux at the wave- 
length of the line. The latter was eval uated for all the 
stars of the spectrophotomctric atlas of iGunn fc Strvker 
f 1983j with the angular diameters calculated through the 
[Barnes fc Evand ( 1976^ relation. The continuum flux of 
HD 171488 was f ound by interpolating the values for the 
Gunn fc Strvked (|l983h stars at B - 1^ = 0.62. As a check, 



we also used the syiithetic low-resolution NextGen spec- 
tra (JHauschildt et al.lll999( ) by interpolating the tabulated 
continuum fluxes at the star temperature. We found nearly 
identical values of the continuum flux (within 5%) with the 
two methods. The average chromospheric line fluxes are 
reported in Table [31 

For evaluating the significance of the ly ™ modulations 
we used simple x^ tests. If the data are totally uncorre- 
lated with the rotational phase, no significant x^ difference 
is expected between a fit of a smooth periodic function and 
that of a constant function (the weighted average of the 



reversal in its core. The Hei A5876 line appears as an 
absorption feature. The Can H&K cores exhibit strong 
emission features and He emission is also visible (Fig. [5])- 

The values of emission equivalent widths, ly™, in the 
different chromospheric diagnostics are listed in Table [2] 
along with their errors. 

We also evaluated the total radiative losses in 
the chromospheric lines following the guidelines by 
Frasca fc Catalanol ( 19941 ). i.e. by multiplying the aver- 



data). As a fitting function, we used a Fourier polyno- 
mial of degree 2 which is able to reproduce asymmetrical 
curves. Additionally, we calculated the probability P^2 ^^2 
that, in a x^ distribution with n — p — 1 degrees of free- 
dom (with n the number of data points and p the num- 
ber of fitting parameters), a random variable is less than 
the measured Xg. A small value of P^2^^2 for the periodic 
function assures a statistically significant modulation. The 
X^ and P^2<^2 values are also reported in Table |31 

5.1. Ha line 

The Ha line is one of the most useful and easily accessi- 
ble indicators of chromospheric emission related to solar 
and stellar activity. Its source function is photoionization- 
dominated in the quiet chromospheres of the Sun and 
solar-like stars, but it can become collision-dominated in 
very active stars or inside ARs. As a consequence. Ha is 
useful for detecting chromospheric solar and stellar plages, 
thanks to their high contrast against the surrounding chro- 
mosphere. 

AU our FOCES spectra of HD 171488 always show 
an Ha absorption profile with a considerable fiUing-in of 
the core. The residual Ha profile is almost symmetric (see 
Fig. [4]) . The values of the residual emission W^ are plot- 
ted as a function of the rotational phase in Fig. [SI where 
a smooth rotational modulation anti-correlated with the 
light curve is apparent. The results of the x^ analysis, 
summarized in Table [21 assure the statistical significance 
of the observed Ha rotational modulation. 



5.2. Hf3 line 

When compared with the non-active template, the H/3 
profile of HD 171488 also appears filled in by core emis- 
sion. This contribution is, however, much smaller than for 
the Ha line, as expected from the lower transition proba- 
bilities for H/3 and the higher photospheric fiux at 4861 A. 
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Fig. 4. Top of each panel: Examples of observed, continuum-normalized spectra of HD 171488 (solid line) in the Ha, 
H/3, HeiDa, and Can IRT regions together with the non-active stellar template (dotted line). Bottom of each panel: 
The difference spectra of the two upper spectra. The hatched area represents the excess emission (absorption for 
HeiDs) that has been integrated to get the net equivalent width. 



The results of the x^ analysis (Table [3]) indicate that 
the small-amplitude variation displayed by W^ 
statistically significant. 



^ is not 



As a further test, we calculated the Spearman's rank- 
correlation coefficient between Wg^ a nd Wg™ by mean s 
of the IDL procedure r_CORRELATE |Press et al.lll986l ). 



We found a very low correlation coefficient p = 6 x 10~ 
with a significance 1.0, which means no correlation. We re- 
mark that this could be due to the faintness of the residual 
emission that prevents us from clearly detecting variations 
over the errors. 

On the basis of the Ha and H/? flux we evaluated a 
Balmer decrement Fug /FnB — 2.7 ± 0.6. As suggested by 
Hall fc RamsevI ( 1992^ . high values of the Balmer decre- 
ment (^ 10) are typical of solar p rominences seen off-limb 



Landman fc Mongillolll979l ). whereas solar plages o r 
flares have always values of about 1-2 (e.g.. lChesteijl991 ). 
Indeed , the theoretical NLTE model developed bv iBuzasi 
(|l989l) leads to high Balmer decrement values, ranging 
from 3 to 15 for prominence-like structures viewed away 
from the stellar disk, while values in the range 1-2 are 
found for plage-like structures. 

The Balmer decrement measured for HD 171488 sug- 
gests that the emission in the Balmer lines is basically due 
to magnetic surface regions analogue to solar plages and 
eventual prominences play a marginal role. 



5.3. Helium D3 line 

The He I A 5876 (D3) line is another useful diagnostics of 
chromospheric activity in the optical domain. It is partic- 
ularly helpful for F-type stars in which the other activ- 
ity indicators in the optical/IR range are more difficult 
to be observed due to t he strong continuum flux (e.g., 
iRachford fc Foighti 120091 ). Both in these stars and in the 
moderately active G an K-type stars , this line is normally 
seen as an absorp tion feature (e.g., iHuenemoerden 1 1 986 : 
iBiazzo et al.ll2007 ). The He i absorption implies a temper- 
ature of ~ 10 000 K in the layer where the line is formed, 
i.e. this line, as in the Sun, is a diagnostics of the upper 
chromosphere. 

The adopted reference spectrum (72 Her) does not 
show any He i absorption nor emission, as expected on the 
basis of its very low activity level. In this case, the spectral 
subtraction technique, cleaning the spectrum from nearby 
photospheric absorption lines, allows to emphasize the He i 
line of HD 171488, which appears as an excess absorption 
in the residual spectra (see Fig. |4]), and to measure its 
equivalent width, WhoI- For W^hcI we adopted the usual 
convention that an absorption line has a positive equiva- 
lent width. 

As seen in Table [5] and Fig. [51 the D3 line is always in 
absorption with small variations only possibly correlated 
with the star rotation. 
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Fig. 5. Top panel: Example of an observed, continuum-normalized spectrum of HD 171488 (solid line) in the Call 
H&K region together with the non-active stellar template (dotted line). Bottom panels: The difference spectrum. He 
emission is barely visible in the raw spectrum, but it is clearly emphasized in the residuals. 



As for the H/?, the x^ test indicates that there is no 
significant rotational modulation of WhcI- However, the 
rank-correlation coefficient with the Ha, p = —0.345 with 
a significance 0.298, indicates a marginal anti-correlation 
between V^g™ and WhoI, in the sense that the Hei ab- 
sorption is slightly weaker when the net Ha emission gets 
stronger. 



5.4. Can H&K lines 

The Can H&K lines show strong emission cores that 
appear nearly symmetric in all spectra without any de- 
tectable trace of self-absorption. He emission is barely seen 
only in the few best spectra and become more evident af- 
ter the subtraction of the non- active template (Fig. [5]). 
Thus, we could not investigate the behavioi ir of He with 
the s tar rotation as we did for SAO 51891 ( Biazzo et al.l 
2009h . 

The absorption wings of each of the two Can lines 
span two eehelle orders and it was necessary to m erge 
them, following the guidelines of iFrasca et al.l ( 20001) . be- 
fore proceeding with the spectral subtraction analysis. The 
non-active template was built up, with the same merging 
procedure, from the FOCES spectrum of 72 Her broad- 
ened at V sin i = 37.1 km s ^ ^ . No emission in the core of the 
Can H and K lines of the template is visible (Fig. [SJ up- 



per panel), as expected from the very low chromospheric 
activity level reported in the literature (see above). 

We measured the equivalent widths by integrating the 
emission profiles in the subtracted spectra, as for Ha. We 
disentangled the contribution of the He line from the Ca n 
H by fitting Voigt functions to both profiles. The sum of 
the net equivalent widths of H and K lines, Wg™j_^^j^, 
displays a significant modulation with the star rotation 
(bottom left panel in Fig [5]) that is not as outstanding 
as that of Ha but its reliability is testified by the results 
of the x^ tests (Table |31). The rank-correlation coefficient, 
p = 0.567 with a significance 0.112, supports a rotational 
modulation of the Can H-l-K emission correlated with the 
Ha one. 

We calculated the radiative chromospheric losses in 
the H&K lines analogously as done for Ha. In particu- 
lar, we used two 10 A- wide bands centered at 3910 and 
4010 A, i.e. at the two sides of the H&K lines, to per- 
form the fiu x calib ration, following the prescriptions by 
Frasca et al.l (J2000r) . The ratio of the residual peak inten- 
sities, tk/^h — 1.08 ± 0.015, is larger than typical values 
measured in the quiet solar chromosphere (rx/^H —1.05) 
but it is smaller t han what it is found in solar plages 
(rK/rn ~1.20, e.g. lLinskvlll970l ). If the solar analogy is 
valid, this would mean that we are observing, at all phases, 
a "mixture" of quiet chromosphere and ARs. 
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5.5. Can I RT lines 

The lines of the Can infrared triplet, which share the up- 
per level of the H and K transitions, present some ad- 
vantages compared to the Can H & K lines. They lie in a 
spectral region with a well-defined continuum, making the 
normalization easier during the data reduction. Moreover, 
they are not significantly affected by telluric lines and are 
less affected by atmospheric extinction than the visible 
and ultraviolet lines. As such, they have become very use- 
ful chromospheric diagnostics in recent years also thanks 
to the high sensitivity of the new detectors to the near 
infrared. 

A filled-in Can IRT-3 (A8862) hne profile is shown by 
iBusa et al.l ( 2007 ) which measured a net equivalent width 
(after the subtraction of a synthetic profile) of 0.511 A 
that is larger than the values measured by us, ranging 
from 0.350 to 0.402 A. The total contribu tion of the three 
IRT lines measured bv lBusaet all (|2007l) in 2002 is 1.46 
A, which is significantly larger than the values measured 
by us (1.01-1.16 A, see Tabled]). This could be due to a 
higher activity level during the Busa et al. observations or, 
more likely, to the different photospheric template used by 
them, that is apparently deeper in its core than our own. 

Measurements of the residual EW of Ca n IRT-3 lines 
(after the subtraction of a solar s pectrum broadened at the 
star's V sin i) are also reported bv ljarvinen et al.l ( 2008[l for 
2001, 2002, and 2005 epochs. They found values compa- 
rable to our determinations in 2001 and 2002 and a larger 
EW, on average, in 2005. Moreover, they detected also a 
modulation of the residual EW anti-correlated with the 
contemporaneous light curve in the last season. 

We found a nice rotational modulation of the net 
equivalent widths of the IRT lines, which becomes more 
evident if we consider the total emission of the triplet 
^Caii-iRT (sse Fig. [nj that closely follows the Ha trend. 
The statistical significance of the rotational modulation is 
assured both from the x^ analysis (Table [3]) and from the 
strong correlation with the W^ curve (p = 0.755 with a 
significance 0.007). 

The fiux ratio ^8542/^8498 = 1-7, indicative of large 
optical depths, is in the range of values (1.5-3) found by 
Chester! (|l99ll ) in solar plages. Solar prominences have 
instead values of ^8542/^8498 ~ 9, typical of an optically- 
thin emission source. This result confirms what we found 
for HD 171488 from the Ha/H;3 flux ratio. 



6. Plage modeling 

The chromospheric line-flux variations contain valuable 
information on the surface inhomogeneities at this atmo- 
spheric level simultaneous to the photospheric diagnostics. 
We have applie d our spot-model IDL code. Macula 



( Frasca et al.ll2005D . to the curves of W^^, W^'^,,_^^r^, and 



wiz 



Ha' "Caii-IRT' 

H+Ki with the aim of further investigating the main 
parameters of the chromospheric plages and the degree of 
spot-plage correlation. Given the scatter in the data, two 



bright spots (plages) are fully sufficient to reproduce the 
observed variations. 

An important parameter that enters in the plage model 
is the flux contrast between plages and surrounding chro- 
mosphere (i^plagc/^chrom)- ValuCS of Fpiago/i^chrom ~ 2, 

that can be deduced averagi ng solar plages in Ha (e.g., 
Ellison! Il952t lLaBonte!ll986bl ). are too low for reproduc- 



ing the ^-^23% variation of W^. In fact, very large plages, 
covering a large fraction of the stellar surface, would be 
required with such a low flux contrast and they do not pro- 
vide a good fit of the observed modulation on the basis of 
the x^- On the other hand, very high values of flux contrast 
(-F^piagc/^chrom ^ 10) would imply very small plages pro- 
ducing top-flattened modulations that are not observed. 
In order to achieve a good fit, a fiux contrast of 5 was 
adopted. 

The solutions essentially provide the longitude of the 
plages, giving only rough estimates of their latitude and 
size. We searched for the best solution starting with the 
values of longitude, latitude, and radii of the photospheric 
active regions and we found that only the longitude of 
the plage at lower latitude must change appreciably. The 
plage radii are, however, strongly dependent on the as- 
sumed flux contrast Fpiagc/i^chrom- This means that only 
the combined information between plage sizes and flux 
contrast, i.e. some kind of plage "luminosity" in units of 
the quiet chromosphere (Lpiagc/^quiot) can be deduced as 
a meaningful parameter. We remark that the ly™ level of 
the quiet chromosphere is also a rather uncertain param- 
eter, because it depends on the plage distribution which 
can give a non-negligible contribution also at the phase of 
minimum W^"^. The solution of the Ha light curve is dis- 
played with a dotted line in Fig. [5] along with the observed 



lyg™ values. 



•^•^Call-H+K 



for the W^^,, 



,IRT and WSZ^ 



tionsoftheT^^™, 



.iRTandI^^-_ 



For the solution of the W^caii-iRT and 
curves we had to use a lower plage contrast with the same 
radii used for the Ha curve, namely Fpiagc/-F'chrom =3.3 
for the former and and 3.1 for the latter. This is consis- 
tent with the variation amplitude of about 14% and 11% 

H+K' respectively. The solu- 
jj_^j^ light curves are also 
displayed with dotted lines in Fig. |6] and the parameters 
of the spots and plages, as recovered by our model, are 
listed in Table U) The values of plage co ntrast in Ha and 
Can are nearly the same found bv Frasca et al.l ( 20001 ) for 
HD 206860, another young {age ~ 300 Myr) early-G type 
star. 

Considering the aforementioned size/contrast degen- 
eracy, this would imply that, in these very active solar- 
type stars, the plages are brighter (compared to the adja- 
cent quiet chromosphere) when observed in the Ha rather 
than in the Can K line. This result is not in contradic- 
tion to what is known from observations of solar ARs in 
UV emission lines, for which the enhancement of inten- 
sity inside plag es increases with t he temperature of line 
formation ('e.g.. iNoves et al.lll970ri . As mentioned above. 



the flux contrast for Ha solar plages, as measured from 
residual line-core intensity, -Fpiagc/^chrom «1.5-2.5 (e.g., 
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Fig. 6. Left panels. From top to bottom. V magnitude, VFg™, W^caii-iRT' ^^'^ ^Caii-H+K versus the rotational phase. 
Note the clear anti-correlation of the photospheric {V light curve) and chromospheric diagnostics. The dotted lines 
represent the solution of spot/plage model. Right panels. From top to bottom. B magnitude, radial velocity, W§[a, and 
Wnei versus the rotational phase. 



Table 4. Spots and plages parameters. 



Diags. 


AR 


Radius 


Lon." 


Lat. 


Pplago 
Fchr 


Tspot 
Tphot 


V 


#1 


15°.5 


10° 


60° 




0.74 




#2 


16°.3 


215° 


75° 




0.74 


Ra 


#1 


19°.2 


50° 


60° 


5.0 






#2 


14°. 


217° 


75° 


5.0 




Ca ll-IRT 


#1 


19°.2 


29° 


60° 


3.3 






#2 


14°. 


217° 


75° 


3.3 




Ca ll-HK 


#1 


19°.2 


29° 


60° 


3.1 






#2 


15°.5 


217° 


75° 


3.1 





'^ Longitude increases with phase, and 0° longitude cor- 
responds to 0.0 phase. 



EUisonl Il952l : iLaBontd Il986bl ) , is much smaller than the 
value derived by us for HD 171488 and HD 206860, while 
pplago /-Fchrom — 3 fouud for the Can p lages falls in the 
range found for the solar plages (e.g., iLaBontd Il986at 
Avres et al.lll986l ) . Anyway, for a proper comparison with 
solar plages, we must consider that for the latter the resid- 
ual photospheric flux in the line core is not removed, and 
indeed it is not easy to be distinguished from the true 
chromospheric emission at the low activity level of the 
quiet Sun. This photospheric contribution is higher in the 
core of Ha than in the H&K Call lines giving rise to a 
contrast drop. To this aim, we made our values of the 
Ha and Ca ii emission equivalent widths comparable with 
the solar determinations by adding back the photospheric 
contribution at the bottom of the lines as measured in the 
non-active template. 



We found a strongly reduced amplitude, ^ 12%, for 
the Ha curve for which our model derives a plage flux 
contrast of 1.8, in line with the average solar values. For 
the Can H&K residual intensity, the modulation drops 
from ^ 11% to ^ 7% which lowers the plage contrast 
from 3.1 to 2.5. This would imply that, for very active 
solar-type stars, the plage contrast (with respect to their 
"quiet" chromosphere) is analogous to what found for the 
Sun but the overall chromospheric flux is much larger and 
the plages are bigger. This gives rise to the wide observed 
line-flux modulations. 

Another interesting result of our plage model is that 
the longitude of the plage at lower latitude (those one giv- 
ing the highest contribution in the flux curves) is larger 
than the one of the corresponding spot by about 20-40°, 
i.e. the plage is preceding the spot in the sense of star's 
rotation. This result is found in all the chromospheric di- 
agnostics for which a significant modulation has been de- 
tected and is confirmed by the simple cross-correlation of 
the light curve with the W^"^ modulations. We found a 
deep minimum of the CCF for the correlation between 
F-light curve and W^ with a phase shift of the Ha with 
respect to the V curve of 0.094 in phase units, which trans- 
lates into a longitude difference of about -1-34°. Slightly 
smaller phase shifts are found for Can IRT (0.070) and 
Can H-l-K (0.052) corresponding to a longitude difference 
of about 20° between plage and spot. 

Lead and lag between plages and spots have been 
suggested to explain phase shifts between photometric 
light curves and chromospheri c ind ic ators that are not 
exactly in anti-phase. Radick et al.l (J1987I ) found spots 
leading plages by about 19° in the Hyades star VB 31 
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(GOV), while ICatalano et al.l ( 20001 ) found both lead and 
lag in the active components of a few RS CVn systems, 
with a tendency for spots to lead plages. The Sun it- 
self, when observed as a star in the integrated light of 
the C II chromospheric line (A 133.5 nm) from UARS 
SOLSTICE experiment, displays alternating positive and 
negative longitude shifts with respect to the spot photo- 



centers (jCatalano et al.lll998l ) 



7. Conclusion 

We performed an accurate analysis of high-resolution 
spectra and contemporaneous BV photometry of the 
nearby young solar-type star HD 171488 with the main 
aim to characterize the behaviour of spots and plages on 
its surface. Revised astrophysical parameters, such as Tcs, 
\ogg, [Fe/H], rotational and heliocentric radial velocity, 
and lithium abundance, were also determined and have 
been found in close agreement with previous determina- 
tions from the literature, with the exception of the metal- 
licity for whi ch we did not find a sub -solar abundance as 
suggested bv lStrassmeier et al.l (|2003[ ). Both the position 
on the HR diagram and the lithium abundance give an 
estimate of about 50Myr for the age of HD 171488, con- 
firming it as a young Sun on its way to the ZAMS. 

We were able to find a clear radial velocity modulation 
with an amplitude of about 500 ms~^ due to a "Rossiter- 
McLaughlin" effect caused by starspots. A spot model- 
ing of the B and V light curves, performed with our code 
Macula, is also able to reproduce the radial velocity mod- 
ulation by adopting only two high-latitude spots 1500 K 
cooler than the quiet photosphere. A more accurate sur- 
face reconstruction was also performed by applying the 
Dopplcr Imaging analysis to two photospheric spectral 
lines and the B and V light curves. The method provides 
maps that essentially confirm the large polar spotted area 
with a temperature difference AT ~ 1500 K with respect 
to the unspotted photosphere. 

Chromospheric activity was also studied with differ- 
ent spectral diagnostics. Equivalent widths of the residual 
emissions in the cores of Ha, Call H&K, and Can IRT 
lines show clear modulations with the rotational phase 
and are ascribed to the presence of active chromospheric 
regions analogue to solar plages. A small but significant 
emission in the H/3 core, as well as an absorption feature 
corresponding to the He i A5876 line, are emphasized by 
the subtraction of the non-active template and are addi- 
tional proofs of the high activity level of HD 171488 also 
in the upper chromosphere. Due to the fairly large relative 
errors, no clear rotational modulation could be detected 
in the H/3 and He i line fluxes. The average B aimer decre- 
ment, measured as the ratio of Ha to H/3 fiux, confirms 
that the observed chromospheric emission is mainly com- 
ing from extended facular regions. 

The Ha and Ca ii fiux curves are well fitted by a model 
with two circular bright spots (plages), also performed 
with Macula, with sizes similar to the underlying spots 
and a flux contrast of 5 for the Ha and about 3 for the 



Can H&K and IRT fluxes, as expected from the differ- 
ent modulation amplitude. The flux curves are nearly an- 
ticorrelated with the photometry, even if the plage at a 
lower latitude precedes in longitude, by about 20°-40°, 
the corresponding photospheric spot. This finding is not 
anomalous since lead and lag between plages and spots 
have been observed in the Sun and are frequently invoked 
as the cause of phase shifts between light curves and chro- 
mospheric line fluxes observed both in young active stars 
and in RS CVn systems. 
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Table 5. BV Johnson photometry of HD 171488. 
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Fig. 7. Observed FOCES spectrum of HD 171488 in the 
Mg lb triplet (upper panel) and 5600 A (lower panel) spec- 
tral regions together with the ELODIE standard spectrum 
broadened at 37.1 km s~^ overplotted with a thin/red line. 
In each box the difference (observed — synthetic) is also 
displayed with a dotted line around zero ordinate. 
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